% M
%
% Files
%   ace                                 - 4-sep-2007. This code reads Julia and ACe data sets
%   ace1                                - 4-sep-2007. This code reads Julia and ACe data sets
%   ace1_lt_dependency                  - 4-sep-2007. This code reads Julia and ACe data sets
%   ace1_phase_rotate                   - 4-sep-2007. This code reads Julia and ACe data sets
%   ace1_prediction_plot                - 
%   ace1_prediction_plot_few            - 
%   ace_amie_joint_analysis             - 
%   ace_burke_test                      - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_esk                             - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_hua                             - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_imf_dependency                  - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_julia_ler                       - 30-nov-2007. This code reads Julia and LER & ESK data sets
%   ace_julia_pcn                       - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_pcn                             - 4-sep-2007. This code reads Julia and ACe data sets
%   ace_pcn_corr                        - 
%   ace_substorm                        - 
%   ace_tf_test                         - load C:\Manoj\projects\ace\20080303\season_all.mat F Txy;
%   add                                 - Add SPMs
%   addspm                              - 
%   ak_cres                             - 
%   amie_julia                          - 4-sep-2007. This code reads Julia and ACe data sets
%   ap2kp1                              - 
%   averagespm                          - code for averaging SPM's
%   avspc                               - function to average spectral estimates
%   AvTrfn                              - function to average transfer functions 
%   bibtex                              - 
%   bin                                 - calculates frequency of occurances
%   bio-nn                              - 
%   BivError                            - function to calculate bivariate errors
%   BPN.M                               - Trains two layer feed forward network with momentum 
%   c_res                               - Compute C-reponse over layered earth
%   cable1                              - Latest series of computation for ocean stuff
%   cable2                              - Copy of cable1.res with modifications to suite the
%   calib05                             - computes the theoretical calbiration
%   Calread                             - clears unneceesary vriables after
%   cchdec                              - Cholesky decomposition of matrix
%   cham_lt                             - calculates the local time from fday and longitude
%   champ                               - read in total magnetic field OHCHAMP. 10.02
%   Champ1                              - read in total magnetic field OHCHAMP. 10.02
%   champ_hua                           - 
%   champ_hua1                          - IMPORTANT NOTE !!!
%   champ_hua_pomme                     - 
%   champ_hua_pomme_rob                 - 
%   champ_hua_temp                      - 
%   champ_hua_temp_rob                  - 
%   champ_hua_test                      - 
%   champ_julia                         - 
%   champ_julia1                        - This file, ie champ_julia1.m was created from champ_julia on 15dec06
%   champ_julia_kp                      - 
%   champ_julia_temp                    - 
%   champ_julia_temp_rob                - 
%   champ_julia_tmp                     - This champ_julia_tmp file was created from champ_julia on 15 dec 06, to
%   champ_lt                            - calculates the local time from fday and longitude
%   change                              - Change block size
%   chavespec                           - Chave's method of making the estimate robust 
%   chebi11_filter                      - Returns a discrete-time filter object.
%   cis                                 - to plot CIS locations
%   cm4_eej                             - generate the variable vector for firing the CM4
%   cmpedi                              - Reads in an edi file
%   cmpedi1                             - Reads in an edi file
%   cmpresph                            - Display app res, phase & asso errors
%   cmprts                              - clears unneceesary vriables after
%   cohe                                - function calculate e-predicted coherency
%   cohguess                            - code for averaging of SPMs wrt coherence
%   cohplot                             - utility to plot
%   cohspm                              - code for averaging of SPMs wrt coherence
%   colorbarf                           - Display color bar for a filled contour plot.
%   compresph                           - Display app res, phase & asso errors
%   compress                            - rts file cmpression
%   CompuSPM                            - Function to compute Cross and auto spectra
%   con2str                             - convert day,hr,sc to string
%   cor_obs                             - Correction for Local Solar Time to IST
%   corpr                               - pearson correlation of 2 ts 			-	15.08.2000
%   CorTre                              - remove trend and bias from time series
%   cpress                              - rts file cmpression
%   createfigure                        - CREATEFIGURE(X1,Y1,Y2,Y3)
%   delflt                              - Delay line filter
%   Despike                             - To remove spikes in time series
%   dispcoh                             - Displays predicted coherencies
%   dispresph                           - Display app res, phase & asso errors
%   disprhc                             - Display app res, phase ,asso errors & coherence, dof et
%   disprts                             - read,display & scroll MT ts data in RTS format. 
%   dispspm                             - Displays APp re phase and errors
%   diw_movie                           - 
%   dmy2jd                              - function to convert 
%   DoNotClose                          - 
%   doparzen                            - multiply parzen window to
%   dummy                               - Process the eej strength derived from satellite data
%   dummy1                              - generate the variable vector for firing the CM4
%   dummy33                             - Written to plot resitivity/phase patches (sections)
%   E_merg                              - FORMAT OF THE SUBSETTED FILE
%   ecco_disp                           - read an display global ocean mag data
%   edi_in                              - read in an EDI file and compute MT app. rs. and phs.
%   edi_out                             - function to write spectral matrices
%   eej                                 - 
%   eej_champ                           - Process the eej strength derived from satellite data
%   eej_champ1                          - Process the eej strength derived from satellite data
%   eej_champ2                          - Process the eej strength derived from satellite data
%   eej_champ3                          - Process the eej strength derived from satellite data
%   eej_champ4                          - Process the eej strength derived from satellite data
%   eej_champ5                          - Process the eej strength derived from satellite data
%   eej_champ6                          - Process the eej strength derived from satellite data
%   eej_champ7                          - Process the eej strength derived from satellite data
%   eej_champ_temp                      - Process the eej strength derived from satellite data
%   eej_cimage                          - function[dx,dz] = eej_dx(w,h,x),
%   eej_climate                         - note 28-feb-08. new climatological model of Patrick was used 
%   eej_dummy                           - Process the eej strength derived from satellite data
%   eej_dx                              - function[dx,dz] = eej_dx(w,h,x),
%   eej_induction_paper                 - Process the eej strength derived from satellite data
%   eej_induction_paper_new             - 9-March-2007 Note: Chnages have been made to the orginal file
%   eej_obs                             - Process the eej strength derived from satellite data
%   eej_period                          - Scrip to do periodociti analysis with Stefan's data
%   eqspec                              - code to compute MT trfn from raw ts files !!!
%   eqts                                - To analyze the time series related to the
%   Esgns                               - Computation of signal and noise power in E fields
%   esig                                - Computation of signal and noise power in E fields
%   esk_hua                             - 4-sep-2007. This code reads Julia and ACe data sets
%   euvac_hua_julia                     - 
%   euvac_hua_julia_rob                 - 
%   fftfrq                              - this code just returns the fft frequencies
%   filt_hua                            - 
%   fischer                             - function to caluculate 
%   florida                             - To compare the transport  estimates from Voltage measurements
%   ga-gm                               - 
%   geeta                               - 
%   gespec1                             - 
%   get_dip_lat                         - double get_dip_lat(double z, double f){
%   get_ut                              - 
%   getbfile                            - simple file open menu
%   getrho1frq                          - 
%   GetRobSpectra                       - function to detrend, de-bias, calibrate, fft, and cross & auto power calculation
%   getrobspectra1                      - function to detrend, de-bias, calibrate, fft, and cross & auto power calculation
%   GetSelArray                         - function to calculate 
%   GetSpectra                          - function to detrend, de-bias, calibrate, fft, and cross & auto power calculation
%   GetTwoSpectra                       - function to detrend, de-bias, calibrate, fft, and cross & auto power calculation
%   GG2GMD                              - Transformation from Geographic Coordinates to Geomagnetic Dipole Coordinates. 
%   globe_plot                          - 
%   grid_min                            - 
%   harita                              - Plot Time series from xls file
%   hua_substorm                        - 
%   hyb_ett                             - Read the HYB and ETT data from 1980 tp present
%   hyd_ett_wdc_correct                 - 18-March-2008. these scripts were used to correct the NGRi wdc files
%   indu                                - 
%   induction                           - ---------Pl do not change the frequency setting !!!!-------------------
%   IsOver                              - function to check overlap of
%   iugg                                - read in total magnetic field OHCHAMP. 10.02
%   JD2DMY.M                            - function to convert 
%   jknf                                - delete one (jack knife) estimate
%   jknfvar                             - delete one (jack knife) estimate
%   jktf                                - Computes bi-variate transfer functions between two channels, given
%   julia                               - 
%   julia_drift_correction_example_plot - 
%   julia_substorm                      - 
%   kak_substorm                        - 
%   Kill50                              - just copy of FN.m
%   kkr                                 - impedance rotation
%   knab                                - To test the 
%   kod_read                            - 
%   kp2ap                               - !caution uses the rounded kp indices (ie .3 and .7 instead .33 and .67)
%   ler_esk                             - 
%   ler_hua                             - 5-dec use esk-hua instead
%   Maps_for_Manoj                      - 
%   medout                              - function to supress outliers
%   medSPM                              - code for averaging SPM's
%   MF_demo                             - 12 AUgust 2005
%   minus1                              - 
%   mtfreq                              - computes target frequency, target lines and radius 
%   MTFW.M                              - MT fwd & inverse computaton - keller & kaufman 
%   mtfwd                               - MT fwd & inverse computaton - keller & kaufman 
%   mtfwd1                              - MT fwd & inverse computaton - keller & kaufman 
%   mtgrid                              - grid  MT profile data stes
%   mtphd                               - stns = [1011,1012,2004,2010,2013,2016,2018,2020,3005,3006,3010,3014,3015,4001,4002,4003,4015,4017,4019];
%   mtproc                              - code to compute MT trfn from raw ts files !!!
%   mtprocm                             - Multiple band processing 
%   mtprocn                             - code to compute MT trfn from raw ts files 
%   mtts                                - Program to read in RTS file, calibrate and write as ascii fil
%   multiplotcolours                    - Specify colors for multiplot data
%   my_patch                            - 
%   myfilter                            - Returns a discrete-time filter object.
%   myfunct                             - 
%   myfunction                          - 
%   mynorm                              - 4-Sep-2007
%   no                                  - 
%   normrnd                             - Random matrices from normal distribution.
%   NormSPM                             - Normalizing SPM wrt power in
%   NrFct                               - --------------------------------------------------------------------------
%   occam_transf                        - This function was written to convert occam rotated coordinates to real
%   oce_mf6                             - write annual harmonics of ocean velocity
%   ocoh                                - computes ordinary coherence
%   parzen                              - Computes Parzen window for current radius
%   pathdef                             - Search path defaults.
%   pchan                               - function to compute EpE
%   pcoh                                - computes predicted coherence
%   period                              - np = number of frequencies examined (= length of output freqs file)
%   phcor                               - phase correlation
%   piu_read                            - 
%   plot_all                            - 
%   plot_ls_regression                  - This function plots regression fit, scatter and returns the correlation
%   plp_champ                           - plp_header = ['1.year|2.month|3.day|4.hh|5.mm|6.ss|7.localtime|8.longitude|9.radius|10.kp|11.smooth|12.uncertainty|...
%   polangle                            - polarization angle (ute)
%   pommecalc                           - This script computes the geomagnetic field values from Pomme
%   PutNoise                            - Adding noise
%   pvmod                               - returns the phase delay at ADC of GMS05
%   qqplot1                             - QQPLOT Display an empirical quantile-quantile plot.
%   r_resid                             - function calculate the residual between
%   r_trfn                              - Computes bi-variate transfer functions between two channels, given
%   rbpz                                - robust parzen averaging
%   rdcomp                              - compression of time series
%   read_cable                          - Reads the *.res files from x3dg (note all the 6 components!)
%   read_char                           - function to read the
%   read_spidr_data                     - 
%   read_tsu                            - File_1 = 'c:\manoj\ocean\x3dg\tsu_results\res2\tsu1.res';
%   readedi                             - Reads in an edi file
%   readedi1                            - read in an EDI file and compute MT app. rs. and phs.
%   readrts                             - code for reading an rts file 18.07.2000
%   resprof                             - reads edi files 
%   rh2ph                               - Estimate the phases from the slope of the apparent resistivity curve
%   rndcheck                            - error checks the argument list for the random number generators.
%   rob_jk                              - code for robust averaging of SPMs
%   rob_regress                         - 
%   rob_regress1                        - 
%   robspm                              - code for robust averaging of SPMs
%   robspm1                             - code for robust averaging of SPMs
%   RrCorrectTS                         - ---------- information ------------------------------
%   rrtransfn                           - Remote reference processig
%   rts_in                              - function for reading an rts file 18.07.2000
%   s-n                                 - computes predicted coherence
%   s_imp                               - MT fwd & inverse computaton - keller & kaufman 
%   s_imp1                              - MT fwd & inverse computaton - keller & kaufman 
%   s_n                                 - Signal power estimate
%   satyam                              - Written for NGRI - Styam Project 'Unified Presentation Framework'
%   scalespm                            - just scaling spectra
%   screen                              - script to test Ute's data screening methods
%   scroll                              - Just to Scroll trhough the SPMatrix (time wise)
%   season_filter                       - Filters fday (in mjd wrt 2000) for seasons
%   SelANN                              - Subroutine to select band 4 (GMS05) time series
%   SetDefault                          - function to set defaults for
%   SetLine                             - Used with tipper calculations
%   settline                            - computes target frequency, target lines and radius 
%   settline_s                          - computes target frequency, target lines and radius 
%   sinu                                - jsinusoids
%   solar_fulx_euvac                    - In the meantime, here is a file of the CHAMP EEJ currents 
%   sortspm                             - 
%   spike                               - stimulate a spike
%   split_wdc                           - 
%   spmap                               - to read in edi files
%   spmap1                              - to read in edi files
%   sq_cor                              - To compute the latitudinal profile of B_x (ionospheric) from CM4 for
%   sq_cor1                             - To compute the latitudinal profile of B_x (ionospheric) from CM4 for
%   sudoku                              - sudoku
%   svs                                 - code to compute MT trfn from raw ts files !!!
%   temp                                - 4-sep-2007. This code reads Julia and ACe data sets
%   temp1                               - 
%   temp2                               - code to compute MT trfn from raw ts files !!!
%   temp3                               - 
%   test                                - 
%   tf                                  - Computes bi-variate transfer functions between two channels, given
%   tf1                                 - Computes bi-variate transfer functions between two channels, given
%   tf2resph                            - Display app res, phase & asso errors
%   tf_dn                               - Down biased (H reference local) transfn 
%   tf_dn_new                           - Down biased (H reference local) transfn 
%   tf_r                                - RR tf
%   tf_up                               - Up biased transfn using E as reference fields
%   tiptrfn                             - Computes bi-variate transfer functions between two channels, given
%   tiwari                              - These scripts were made to get
%   transfn                             - Computes bi-variate transfer functions between two channels, given
%   tsu_pre                             - To interoplate and fft the tsunami v elocities
%   tsu_res                             - To read and show the results of TS simulation
%   tsu_res1                            - To read and show the results of TS simulation
%   tsu_wav                             - 
%   usefulscripts                       - Hei some nice scripts for auto open, download urls
%   ute                                 - script to test Ute's data screening methods
%   wdc_read                            - read the wdc hourly means
%   WeightSPM                           - returns the weighted sum of spectra
%   workshop                            - 
%   wr_ann_oce                          - write annual harmonics of ocean velocity
%   wr_ann_oce_test                     - write annual harmonics of ocean velocity
%   wr_cable                            - To write the source current files 
%   wr_cable1                           - To write the source current files 
%   wr_cable_temp                       - To write the source current files 
%   wr_oce1                             - To write the source current files 
%   wr_polar                            - To write the source current files 
%   wr_tsu                              - To write the source current files 
%   wr_tsu1                             - To write the source current files 
%   wr_tsu2                             - Used to make source file with Neetu's new data
%   yap_cor                             - 
%   yap_cor_new                         - YAP data for 2002 June 9-14 available at
%   yap_cor_temp                        - 
%   BPN.M                               - 
%   BPN.M                               - 
%   JD2DMY.M                            - 
%   BPN.M                               - 
%   JD2DMY.M                            - 
%   MTFW.M                              - 
%   bio- 
%   ga- Constrained optimization using genetic algorithm.
%   s- 
%   MTFW.M                              - 
